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THE LARI METHOD FOR ISO-CAM/PHOT DATA REDUCTION AND ANALYSIS 
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ABSTRACT 

The techniques and software tools developed for the 
reduction and analysis of ISO-CAM/PHOT data with the 
LARI method are presented. The method, designed for 
the detection of faint sources in ISO raster observations, 
is based on the assumption of the existence of two dif- 
ferent time scales in the detectors' transient behaviour, 
accounting either for fast or slow detectors' response. 

The specifically developed IDL software includes: a 
reduction pipeline performing basic operations such as 
deghtching and background determination; the fitting pro- 
cedures proper, modelling the time history of individual 
pixels and detecting any fiux excess with respect to the lo- 
cal background ascribable to potential sources; mapping, 
source extraction and fiux estimation procedures; simula- 
tion procedures allowing one to estimate the errors arising 
from different instrumental and reduction effects. More- 
over, an easy-to-use graphical user interface allows one 
to quickly browse the data and carry out the substan- 
tial amount of interactive analysis required when the au- 
tomatic fit fails and to check the reliability of detected 
sources. 

This method provides source lists of great reliability 
and completeness and an outstanding photometric accu- 
racy, particularly at low redundancy levels, where the re- 
liability of ISO-CAM/PHOT source lists, even at moder- 
ately bright fiux levels, has been a long-standing issue. 

In this work a description of the techniques and soft- 
ware tools that have been developed is given, alongside 
with some highlights from the results obtained thanks to 
their application to different fields. 
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1. Introduction 

All data gathered by the ISO satellite, and particularly 
those from the its two cameras, ISO-CAM and ISO-PHOT 
(hereafter simply CAM and PHOT), are very difficult to 
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reduce, both due to the strong transient behaviour of the 
cryogenically cooled detectors, e .g. after a change in the 
incident fiux ( jCoulais et al. 200C ) , and to the frequent and 
severe cosmic ray impacts yielding qualitatively different 
effects (known a s common glitches , faders, dippers, drop- 
outs and others, Claret et al. 1998 ). 



While it was variously demonstrated that it is possi- 
ble, at least to a certain extent, to satisfactorily describe 
the satellite's different detectors' behaviour adopting some 
physical model, the large number of readouts involved in 
raster observations and the peculiar nature and strength of 
noise patterns also require efficient and robust algorithms 
to be developed so as to make the actual data reduction 
undertaking feasible in a nearly-automatic way. 

A number of data reduction methods has thus been de- 
veloped and tested, mostly on CAM deep fields (e.g. the 
PRETI method bylStarck et al. 19991 and the Triple Beam 



Switch method by Desert et al. 199£). Unfortunately, such 
methods proved useless for all PHOT data or on CAM 
shallower fields, leading to a high number of false detec- 
tions and severe incompleteness. Besides, these methods 
suffered from the lack of an efficient way to interactively 
check the quality of the data reduction when needed. 

The LARI method, originally developed as an answer 
to the problems posed by the reduction of the obscrvatons 
carried out as part of the ELAIS survey (Oliver et al. 2000) 
and first presented in Lari et al. 2001, has been devised to 



overcome these difficulties and provide a fully-interactive 
technique for the reduction and analysis of CAM/PHOT 
raster observations at all flux levels, particularly suited 
for the detection of faint sources and thus for the full ex- 
ploitation of the scientific potential of the ISO archive. 

2. The Model 

The LARI method describes the sequence of readouts, or 
time history, of each pixel of CAM/PHOT detectors in 
terms of a mathematical model for the charge release to- 
wards the contacts. Such a model is based on the assump- 
tion of the existence, in each pixel, of two charge reser- 
voirs, a short-lived one Qb (breve) and a long-lived one 
Qi (lunga), evolving independently with a different time 
constant and fed by both the photon flux and the cosmic 
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rays. Such a model is fully conservative, and thus the ob- 
served signal S is related to the incident photon flux / and 
to the accumulated charges Qb and Qi by the 



g _ J _ ^Qtot J dQb dQi 



(1) 



dt dt dt 

where the evolution of these two quantities is governed 
by the same differential equation, albeit with a different 
efficiency and time constant 
dQ^ 



dt 
so that 



Cil — ai where 



b,l , 



(2) 



S ^{l-eb-ei)I + abQl + aiQf . (3) 

The values of the parameters and ai are estimated from 
the data and are constant for a given detector, apart from 
the scaling of a.i for the exposure time and the average 
signal level along the pixel time history, which is governed 
by the 

= r V I" "^'0 ' ('^) 

to V OQ 

where ai^ is the value of ai relative to a reference exposure 
time to and average signal level Sq. The model for the 
charge release, however, is exactly the same for CAM and 
PHOT detectors. 

In practice, an additive offset signal attributable to 
thermal dark current (which is not accounted for in CIA/PI A 
pipeline dark current subtraction) , is added to both S and 
/ in the equation above when it is estimated to be impor- 
tant, i.e. when the deepest dippers' depth exceeds 10% of 
the background level. 

The glitches (i.e. the effects of cosmic ray impacts on 
time history) are identified through filtering of the time 
history and modelled as discontinuities in the charge re- 
lease, leaving as free parameters the charges at the begin- 
ning of the time history and at the peaks of glitches. 

Iteration of the fitting procedure is interrupted when 
either a satisfactory data-model rms deviation is achieved 
or the maximum number of allowed iterations is reached. 
In Figure |l] it is shown how a successful fit is thus able to 
recover useful information (specifically, source fluxes) from 
otherwise troublesome parts of the pixel time history. 

3. The Method 
The reduction pipeline consists of the following steps: 



PHOT ramps' linearization (following Rodighiero el 
al. 2001) and data smoothing (i.e. median averaging 



over a suitable number of readouts) 

— Standard CIA/PI A raster structure and data-reduction- 
dedicated liscio structure building 

— Dark current subtraction, stabilization (or global, as 
opposed to local, see below) background level estima- 
tion, bright sources' and glitches' identification 

— Time history fitting procedure and interactive " repair" 
on fitting failures 



— Interactive checks on sources detected in time history 

— Flat-fielding, mapping (i.e. projection of fluxes mea- 
sured along the pixel time histories onto a sky map), 
and source extraction (using DAOPHOT's find, par- 
ticularly suited for the detection of point-like sources, 
as implemented in IDL Astronomy User's Library) 

— Interactive checks on sources detected on sky maps and 
back-projected on pixel time histories 

— Source flux autosimulation 

The fitting procedure describes the time history of indi- 
vidual pixels according to the mathematical model seen 
in Section ||, allowing the determination of the breve and 
lunga charge levels, the local background (i.e. the signal 
to be expected on the basis of the previously accumulated 
charges if only the stabilization background flux were hit- 
ting the detector) and the flux excess ascribable to poten- 
tial sources at any given readout). 

The so-called autosimulation procedure for source flux 
estimation accounts for mapping effects in the determi- 
nation of the total flux of detected sources through the 
following steps: 

— First guess of source total flux, based on its observed 
peak flux on the map 

— Back-projection of source from sky map on pixel time 
histories, i.e. construction of theoretical pixel time his- 
tories 

— Projection of theoretical pixel time histories onto a sky 
map, i.e. construction of a theoretical sky map 

— Determination of theoretical peak flux on theroretical 
sky map 

— Source flrts-guess total flux correction based on ob- 
served / theoretical peak flux ratio 

Other factors affecting the source flux estimates, namely 
those arising from the detectors' transient behaviour and 
possible systematic deviations from nominal sensitivities 
as well as from the reduction technique, are then evaluated 
through simulations and absolute flux calibration. 

Once the reduction of all rasters of interest has been 
completed according to the recipe above, one can deter- 
mine the necessary corrections to nominal astrometry (e.g. 
performing a cross-correlation analysis between the list of 
detected sources and a suitable external catalogue) and 
project nearby or repeated fields onto a common mosaic 
sky map, on which source extraction, autosimulation and 
interactive checks can furtherly be performed so as to in- 
crease the quality of the reduction through cross-checks of 
detected sources on different rasters, thus partly overcom- 
ing the severe problems at their boundaries. 

4. The Software 



The method rehes on CIA (|Ott et al. 2001] ) and PIA 
(Gabriel and Acosta-Pulido 199!; ) for basic raw data read- 
ing and manipulation and on home-made IDL routines for 
the data reduction proper. The massive necessary amount 
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of interactive analysis is carried out with an easy-to-use 
graphical user interface, a screenshot of which is shown 
in Figure ^, which allows any kind of visual check and 
interactive "repair" which may be necessary. 

5. Results / Work in Progress 

All parameters indicating the goodness of the reduction 
(reliability, completeness, astrometric and photometric ac- 
curacy) are heavily dependent on the adopted observing 
parameters as well as on the thresholds chosen in the in- 
teractive "repair" , and thus can only be evaluated through 
simulations. 

While the catalogue resulting from the reduction of 
the first portion of ELAIS 15 fim data (Lari et al. 2001) 
is successfully being used for different purposes (see e.g. 



Gruppioni et al. 2002 and Matute et al. 2002), and the 



reduction of several different fields has already been com- 
pleted, simulations and accurate photometric calibration 
are still being carried out, so that it is not presently possi- 
ble to show detailed results. A list of the different projects 
being carried out includes: 

- ELAIS 15 Atm and 90 fim fields ( [Vaccari et al. 2002D 

— Lockman Hole Shallow (LHS) and Deep (LHD) 15 fim 



and 90 fiui fields (Fadda et al. 2002 and Rodighiero et 
al. 2002| , see also Figure g) 



— Hubble Deep Field North and South (HDFs) 7 fim and 
15 ^ui fields 

— A few nearby galaxy cluster 7 /im and 15 fim fields 

while highlights from the expected results can thus be 
summarized: 

— A catalogue of around 2000 15 sources in the 0.3- 
100 niJy flux range from LHS and ELAIS fields 

— Catalogues of the uttermost quality in smaller, cosmo- 
logically relevant fields such as the LHD and the HDFs 

— Flux-level-dependent photometric calibration based on 
predicted stellar fluxes (CAM) or on internal/external 
calibrators' reduction (PHOT) 

— Unambiguous comparison of LARI fluxes with those 
obtained with different methods, e.g. on HDFs 



6. Conclusions 

Originated as an answer to the problems posed by the 
ELAIS data reduction, the LARI method has evolved into 
a complete and well-tested system for ISO-CAM/PHOT 
data reduction and analysis, especially designed for the 
detection of faint sources and the interactive check of de- 
tected sources. Raster observations carried out with ISO- 
CAM LW detector at 7 and 15 /im and with ISO-PHOT 
ClOO detector at 90 fim have been successfully reduced, 
while tests are foreseen to apply the method to other ISO 
detectors, to ISO-PHOT C200 in particular. 

Interactive by its very nature, the method both allows 
ISO-CAM/PHOT data reduction at all flux levels from 



scratch and to check the quality of any independent data 
reduction undertaking, thus leading to extremely reliable 
and complete source catalogues. It is thus believed that 
the LARI method can prove a very efficient tool in pro- 
viding the community with long-awaited agreed-upon re- 
sults from ISO extragalactic surveys, possibly the ultimate 
legacy of the ISO mission. 
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Figure 3. 15' x 15' Lockman Hole 15 fim map, with over- 
Figure 1. Different troublesome situations m CAM pixel ^^^^ gg contours (blue lines) and radio sources (green 
time histories: a) Recovery of stabilization background i^^^^gi^s^ f^^tm de Ruiter et al. 1997) 
level after a fader b ) Recovery of stabilization background 
level after a dipper c) Detection of a bright source hid- 
den by a strong common glitch d) Faint source hidden by 
the recovery of the stabilization background level after a 
dipper. 



